Critical Survey of Isoscalar and Isovector Contributions to the Spin Orbit Potential in 

Relativistic Mean Field Theory 
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The spin-orbit (SO) interaction, emerging naturally from the Relativistic Mean Field (RMF) 
theory is examined critically in the light of the recently measured excitation energy differences 
between the terminating states built on two different configurations for nuclei belonging to the 
lower pf shell. The calculations are carried out using the cranked RMF framework. To further 
probe the iso- vector dependence of the spin-orbit potential, the energy spacing between the gr/2 
and hn/2 states in the Sb-chain is compared to experiment, ft is found that the calculation at the 
quantitative level deviates strongly from the experiment, fn particular the balance of the iso-scalar 
and iso-vector strengths of the effective one body SO potential indicates that additional terms like 
tensor couplings may be needed to account for the experimental data. 

PACS numbers: 



Low energy excitations and global properties in atomic 
nuclei are rather well described by the modern effective 
forces. Different realizations of the short range nuclear 
interaction have been developed, like the Skyrme and 
Gogny for non relativistic Hartree Fock potentials and 
the Relativistic Mean Field (RMF) Lagrangian, based on 
meson couplings. These effective theories employ a set of 
parameters that are adjusted to selected experimental ob- 
servables in order to account for basic nuclear structure 
properties. Since there exist a multitude of observables 
and no real consensus about which of these are to be 
selected in a unique manner, there exist a multitude of 
different force parameterizations. One of the basic prob- 
lems in adjusting the force parameters is related to the 
fact that the single particle states, that are so crucial 
for high accuracy calculations, are in general coupled to 
collective motion and therefore difficult to determine. In 
this context, any dataset representing single-particle mo- 
tion manifests itself as an invaluable source for a rigorous 
test and subsequent fine-tuning of the parameters. 

Recently, it was shown 0, Q that a set of terminating 
states in the mass region A~ 45 may provide a unique 
tool to determine the strength of the spin-orbit (SO) in- 
teraction. The long experience in the studies of terminat- 
ing states within the Nilsson model have clearly shown 
that these belong to the purest single particle states avail- 
able in nuclear structure 0, 0] and hence provide an ex- 
cellent laboratory for the adjustment of effective forces. 
Indeed, in Refs. [l], 0] it was demonstrated that for some 
of the Skyrme force parameterizations, the SO interac- 
tion can be adjusted in such a way, that it agrees with the 
experimental data, whereas some other force parameter- 
izations can be ruled out as irreparable. The robustness 
of terminating states in determining properties of the nu- 
clear energy density functional was further supported by 
the recent comparative study between the fully corre- 



lated shell model and the Skyrme Hartree-Fock (SHF) 
|5j, showing essentially a one to one correspondence. 

Since the SO potential emerges naturally within the 
RMF Q , it becomes highly interesting and important to 
carry out a critical survey of the SO interaction within 
the RMF framework, in the light of the recently measured 
terminating states in the lower pf shell. 

In addition, there is a long standing debate on the iso- 
vector strength of the SO potential for different Skyrme 
parameterizations, exemplified in the recent extension Q 
of the Skyrme SO potential, which needs to be confronted 
to experiment. Similarly, in the RMF, there is an open 
question as to what extend the Hartree potential that 
generates the SO field is sufficient to account for the 
experimental data or whether additional terms or con- 
straints, e.g. tensor coupling [8[ need to be added in the 
theoretical framework. 

One may also consider the Fock terms involving the it 
meson, which strongly affect the effective one body SO 
potential In this work, we carry out a systematic 
investigation of the SO potential by comparison to the 
recent data from the lower pf shell within the framework 
of the cranked RMF (c-RMF) [HOD]- To further address 
the iso- vector dependence of the SO potential, we employ 
the recently studied Ch/2 ~ ^11/2 spacing, E gh , in the 
chain of Sb isotopes 

Recently, experimental data on the terminating states 
in the lower pf shell has been reported in the literature 
03 OS 03 03 03 HI- The terminating states 



considered here involve the fy 2 as well as the d 3 ^ 2 
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configurations where n denotes the number of valence 
particles outside the 40 Ca core and the particle hole ex- 
citation across the magic gap 20 is of proton type. Sum- 
mary of the relevant experimental results is presented 
in Table |TJ Following the discussion in Refs. 0, [H, 
we extract the difference, AE, between the excitation 
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TABLE I: Summary of experimental results on terminating states in lower pf shell. 



Nucleus 


J7/2 


,-1 ,n + l 
"3/2-' 7/2 


AS (MeV) 


Imax 


E max (MeV) Ref. 


Imax 


S ma:c (MeV) Ref. 


42 Ca 


6+ 


3.189 


ri3] 


IV 


8.297 


[13] 


5.108 


44 Ca 


8+ 


5.088 


[14] 


13~ 


10.568 


M 


5.480 


44 Sc 


11+ 


3.567 


[15] 


15" 


9.141 


[15] 


5.574 


45 Sc 


23/2" 


5.417 


[16] 


31/2+ 


11.022 


[16] 


5.605 


45 Ti 


27/2" 


7.144 


[17] 


33/2+ 


13.028 


[18] 


5.884 


46 Ti 


14+ 


10.040 


[19] 


17" 


15.550 


[19] 


5.510 


47y 


31/2" 


10.004 


[20] 


35/2+ 


15.259 


[20] 


5.255 



energies of the states terminating within the d^^f™^ 1 
and fy 2 configurations. The value of AE depends pre- 
dominantly on the size of the magic gap 20 which in 
turn directly relates to the strength of the SO potential. 
Within the spherical Nilsson Hamiltonian [2l| the magic 
gap Ae2o = hhj (l — 6k — 2np), depends on three major 
factors: i) the bulk properties of the potential character- 
ized by huj , ii) the strength of the SO term n and iii) 
the flat-bottom and surface properties entering through 
the orbit-orbit term, ~ p. All three terms influence the 
magic gap 20 with a well defined hierarchy. For light 
nuclei one can disregard the flat-bottom effect (p ~ 0). 
Note that the smaller the energy spacing between 
and I/7/2, the larger the effective SO potential. 

The width of the Nilsson potential, hoj Q , is fitted to nu- 
clear radii. Hence, it sets a common energy scale for all 
microscopic models in low energy nuclear physics. Con- 
sequently, any adjustments of, in particular, the RMF 
parameters leading effectively even to small variations 
in huj in light nuclei is expected to endanger the good 
agreement already achieved between the theory and the 
experiment for radii and binding energies in heavy nu- 
clei. Therefore, the difference in the excitation energies 
provides a unique testing ground for the SO interaction 
usually employed in the mean field calculations. 

These simple arguments are strongly supported by the 
self-consistent SHF calculations [J, [§]■ In particular, 
in Ref. [l[ it is shown that the difference between the 
SHF and the corresponding experimental values of AE, 
averaged over all the available data, i.e. the quantity 
which measures the quality of the parametrization is di- 
rectly correlated with the isoscalar-effective-mass scaled 
strength (true strength) of the iso-scalar SO interaction. 
This result not only exemplifies the role of the SO po- 
tential but also shows that for physical observables like 
AE (or E 9h ) the impact of non-local effects (effective 
mass) on the HF single-particle spectra is, as expected, 
compensated by the effective coupling constants in the 
energy density functional. 

In the RMF theory, the point nucleons are assumed 
to be interacting only via the c-number electromagnetic 
(e.m.) and the a, u> and p meson fields. Within the 
standard nonlinear (a, u>, p) interaction Lagrangian [22I 
[23j the SO potentials for protons and neutrons (U[ s and 



UJl) can be expressed in terms of iso-scalar and iso- vector 
contributions as (see, for example, [24T|): 



[(Cl + Cl)V rP {r) 
±C 2 p \7 rPpn (r)] 



(1) 



Here, m (m*) is the nucleon (effective) mass; the con- 
stants Ci (i = a,uj,p) are defined as C, — mgi/rrii, mi 
being the meson masses; p(f) is the total nucleon density 
(sum of neutron and proton densities) and p pn {f) is the 
difference between proton and neutron densities. Sign + 
(-) arises for protons (neutrons). Clearly, the iso-scalar 
term of the SO interaction is dominant in comparison 
with the iso-vector term. This can be substantiated fur- 
ther by examining the ratio of iso-scalar to iso-vector 
strengths of the SO potential in Equation JT]). For NL1 
[25| and NL3 [26] Lagrangian parameter sets, this ra- 
tio turns out to be 16.54 and 19.62, respectively. For 
Lagrangian parameter sets with density dependent cou- 
plings, DD-ME1 H3 and DD-ME2 [Hj], this ratio is typ- 
ically around 20 or even larger. Further insight into the 
properties of the SO potential can be gained by assuming 
that proton and neutron densities have the same shape. 
In this case, Equation (JTJ) reduces to: 



p(n) _ -{Cl + CpA 



rm 2 m* 



it a; 



N — Z 



A 



9{r) (2) 



Here, — g(r) is the derivative of the density, and is neg- 
ative definite; and A = C p /(C 2 + C 2 ). As expected, 
the spin orbit potential decreases (increases) for protons 
(neutrons) as N — Z increases. However, in contrast to 
non-relativistic models, the iso-vector strength, A, is at 
least an order of magnitude smaller, and is not an inde- 
pendent quantity that can be adjusted further. 

In this work, we solve the c-RMF equations using the 
basis expansion technique (22l |. The upper and lower 
components of Dirac spinor as well as the meson and e.m. 
fields are expanded in terms of three dimensional har- 
monic oscillator basis. We use 8 fermionic and 10 bosonic 
shells. We have explicitly checked for a few representative 
cases that the use of 10 fermionic and 10 bosonic shells 
produces essentially the same values of AE. Since pri- 
mary focus is on relative quantities like AE, inclusion of 
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a larger number of shells does not alter the conclusions 
drawn. Further, pairing correlations are ignored. This 
is justified, since it is expected that the pairing correla- 
tions are quenched for the states considered here [1, 0] • 
Currents are crucial for the calculation of AE and are 
taken into account appropriately. For further technical 
details of the c-RMF model, see JToL [Tl| . The calculation 
are carried out usin g tw o frequently used parameter sets, 
NL1 [H and NL3 [26[. In the calculation of E gh , we 
also employ the RMF model with density dependent me- 
son - nucleon couplings l27l.[28l . l29j using the Lagrangian 
parameter set DD-ME2 28] . Preliminary results of the 
present calculation may be found in Ref . [30l ] . 
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FIG. 1: The difference between theoretical and the corre- 
sponding experimental values of the excitation energy differ- 
ences (AE 1 ) for the terminating states within the d^, 2 f"X 
and fy 2 configurations. 

The difference between the calculated and the corre- 
sponding experimental values of AE is plotted in Figure 
[TJ It is clear that AE values obtained using the param- 
eter sets, NL1 and NL3, though qualitatively in agree- 
ment with the experiment, differ from each other and 
from the experiment. Both over-estimate the experiment. 
The deviation is in general larger than the corresponding 
Skyrme calculations. Both NL1 and NL3 values exhibit 
similar systematics and differ from each other roughly 
by 0.5 MeV on the average. Further, the NL1 values 
are found to be closer to the experiment than those ob- 
tained by using set NL3, with mean squared deviations 
0.59 and 1.31 respectively. Apparently, the magnitude 
of the iso-scalar spin orbit potential is at variance with 
experiment. It should be noted that the nuclei in the 
above figure are arranged in the sequence of increasing 
(N-Z)/A. Hence, a slope in AE is related to inadequacy 
of the iso- vector strength. An inspection of Figure [1] re- 
veals that for both the parameter sets, the discrepancy 
between experiment and theory is maximum for 47 V and 



minimum for 44 Ca. Thus, the discrepancy decreases with 
(N-Z)/A. This analysis points towards deficiency of the 
SO interaction, especially, in two ways: i) the iso-scalar 
term is obviously too weak ii) the iso- vector term requires 
a substantial increase. 

To further investigate the iso-vector part of the SO 
potential, we study the chain of Sb isotopes in the RMF 
framework. It should be noted that for this particular 
case, we solve the RMF equations without taking the 
currents into account. It has been shown [3l| that the 
effect of time reversal breaking which exists in the case 
of odd A systems seems to cancel while taking the dif- 
ferences between single particle energies. Therefore, the 
conclusions drawn here will not change, even if the cur- 
rents are explicitly taken into account. 

The Sb chain gives an important clue to the iso-vector 
dependence of the SO potential from the spacing between 
single proton 2g 7 / 2 and l/in/2 states, revealing that the 
spacing E gh increases with neutron number, i.e the ef- 
fective SO potential for protons decreases with neutron 
number. Note that the single particle states 2g 7 / 2 an d 
lh n / 2 are less pure than the terminating states discussed 
above. For example, the energy spacing between g 7 / 2 and 
h\\/2 states is also influenced by quadrupole collectivity, 
particularly in the mid - shell region, yielding solutions 
that are somewhat deformed. Therefore, these states do 
not allow an adjustment at the same level of accuracy. 
Still, following the discussion in Ref. [12|, the spectro- 
scopic factors show that these states are of dominant 
single particle character, implying that they indeed can 
serve to probe both iso-scalar and iso-vector strength. 

One should mention here that the increase of E gh 
versus N is also discussed in the literature within the 
context of the nuclear shell-model, as being caused by 
the monopole part of the proton-neutron interaction and 
is commonly dubbed as monopole migration effect, see 
e.g. [12] .According to recent investigations, the source of 
monopole migration is a two-body shell-model tensor in- 
teraction [33| . The tensor interaction can naturally be in- 
corporated into the energy density functional framework 
[l,[34j]- In fact, the tensor interaction has been included 
in the tensor interaction in the recent Skyrme - Hartree 
- Fock studies [H, [H, [37j ■ It has been shown there that 
it plays an important role in the evolution of the shell 
structure, for example, the observed trend in the 2g7/2 - 
l/ill/2 spacing. It is well known that in the mean-field 
approximation the two-body tensor interaction leads to 
a one-body SO potential, i.e. it works effectively like the 
two-body SO interaction. The basic difference between 
these two sources of one-body SO potential is that the 
form factor generated by the two-body tensor interaction 
depends on the shell-filling unlike its counterpart com- 
ing from the two-body SO interaction. In an analogous 
manner, the pNN tensor coupling has been discussed 
for the RMF [8| generating an additional contribution to 
the effective one body SO potential. Hence, within the 
mean-field approximation, the monopole migration can 
be interpreted as a dynamical manifestation of an effec- 
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tive one-body SO potential. For a pedagogical discussion 
on this issue, see Ref. [38| . Consequently, with a suitable 
parameterization of the SO potential, one should be able 
to reproduce empirical values of E gh . In Figures [2] and [3] 
we demonstrate that this is indeed the case by showing 
the results of a microscopic-macroscopic calculation with 
a Woods-Saxon potential, including an iso-vector term 
with the same sign as in RMF but considerably larger 
relative strength (see, for example, [39^. 




110 115 120 125 130 135 
Mass Number 

FIG. 2: The calculated and the corresponding experimental 
[l2| E 9h values for Sb isotopes. Solid line marks the result 
of microscopic-macroscopic calculations with a Woods-Saxon 
potential. 

The calculated and the corresponding experimental 
[l2| values of E gh are presented in Figure [2] for all the 
Lagrangian parameter sets considered here. The figure 
reveals that all the calculations deviate by large amounts 
from the experiment. The experimental E gh varies be- 
tween 0.5 MeV to around 3 MeV, whereas the corre- 
sponding calculated values (NL1 and NL3) vary within 
1 MeV only. DD-ME2, on the other hand, yields almost 
a constant spacing. Note that the larger the calculated 
splitting, the smaller is the effective SO potential. The 
difference between the theoretical and the corresponding 
experimental values of E gh is presented in Figure [3J This 
figure highlights the large deviation between theory and 
experiment. Interestingly, it is found that the difference 
decreases with increasing mass number (hence with in- 
creasing (N-Z)/A), for all the three parameter sets. This 
is similar to the case of AE (see Figure [I} . 

We notice again, that the effective iso-scalar SO po- 
tential in the RMF theory is too weak, yielding a much 
too large spacing between the 2g 7 / 2 and I/111/2 states. 
Since the iso-vector dependence of the SO is very weak 
in the RMF theory, with increasing (N-Z) /A, the experi- 
mental values approach the calculated ones. By increas- 
ing the coupling constants of the a and u> meson, one 
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FIG. 3: The difference between the calculated and the corre- 
sponding experimental [l2l ] E gh values for Sb isotopes. 



may increase the effective SO potential, and at the same 
time keep the central potential unchanged. However, this 
will only reduce the relative iso-vector strength of the 
SO interaction further, see Equation (|2|). In this way, 
one indeed is able to improve the effective SO potential, 
however, at the same time, the iso-vector potential be- 
comes weaker. Our earlier study of the symmetry energy 
[4fj| revealed that the coupling strength of the p meson 
is rather well adjusted to reproduce the strength of the 
iso-vector central field. Hence, our results suggest that 
the presently used standard RMF Lagrangians do not 
allow to simultaneously adjust the iso-vector strength, 
i.e. increase its value in agreement to the observed ex- 
perimental values and at the same time increase the SO 
potential. Recently, Piekarewicz [4l[ has studied the sin- 
gle particle spectra near A = 40 by adjusting m a and g a 
to reproduce the ld 3 / 2 ~ ^ s i/2 energy gap in 40 Ca within 
the conventional RMF approach. It should be however 
mentioned that this is a purely localized study around A 
= 40. With only a few parameters in the conventional 
RMF, it seems rather difficult to get a good global fit 
that yields the correct single particle spectrum, like the 
gap lcZ 3/ 2 - 2s 1/2 or spacing lh 11/2 - 2g 7/2 . 

To summarize, a systematic study of the energy split- 
ting of the terminating states between fy 2 and d^^f™^ 1 
configurations using the cranked Relativistic Mean Field 
theory is carried out. The principal aim of this work 
is to critically examine the SO potential, which emerges 
naturally from the RMF framework. Calculations have 
been carried out for NL1 and NL3 parameter sets. It is 
found that the one-body SO potential is essentially too 
weak and at the same time the experimentally observed 
iso-vector dependence stronger than what is found in the- 
ory. The analysis of the energy spacing between the 2g 7 / 2 
and l/in/2 states in Sb isotopes supports these findings. 
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The iso-scalar part of the SO interaction is so dominant, 
that there is very limited freedom to improve the iso- 
vector SO potential generated by the p meson. Even the 
inclusion of density dependence of the meson - nucleon 
coupling strengths does not seem to cure this problem. 
Thus, it can be concluded that the standard version of 
the Lagrangian used, does not give enough freedom to 
achieve the desired strength of the SO potential and the 
correct iso- vector dependence in the SO interaction. This 



may point to the need to extend the RMF model, for ex- 
ample, by including tensor couplings Q or to extend the 
scheme towards Hartree-Fock, including 7r-mesons |9J. 
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